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Abstract 

We consider an inverse problem of identifying the unknown cavities in a heat con- 
ductor. Using the Neumann-to-Dirichlet map as an input data, we develop a linear 
sampling type method for the heat equation. A new feature is that there is a freedom 
to choose the time variable, which suggests that we have more data than the linear 
sampling methods for the inverse boundary value problem associated with EIT and 
inverse scattering problem with near field data. 

1 Introduction 

The biggest advantage of the linear sampling method for inverse scattering problem is that 
it is a very easy numerical scheme to implement. It has attracted and still attracting wide 
audiences in inverse scattering problem. Nevertheless, the validation of the scheme has 
some problem and it was compensated to some extent by the works of Arens et al. ([T], 
[2]). The linear sampling method for inverse boundary value problem for EIT (electric 
impedance tomography) was first given in the unpublished paper by Somersalo f|18j) to 
identify unknown inclusions. A similar linear sampling method was considered by Cakoni 
et al. (0]) for inverse scattering problems using many near field measurements. In this 
paper, we will develop a linear sampling type method for the heat equation which could be 
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of some interest if we not only observe something similar as in the case of inverse boundary 
value problem and inverse scattering problem, but also something new. We will also call 
this method the linear sampling method. 

Of course we cannot have scattering problem for the heat equation. A natural inverse 
problem for the heat equation would be an inverse boundary value problem which conducts 
measurements at the boundary and identify any unknown objects such as inclusions, cavities 
and cracks. In this paper we will establish a linear sampling method to identify unknown 
cavities inside a heat conductor by many boundary measurements. 

There are many other schemes to identify the unknown objects in heat conductors. We 
only give some of the references for these schemes. They are the dynamical probe method 
by Nakamura et al. ([7], [5], [T3], [IS], [12]) and enclosure method by Ikehata et al. ([5], 

m, m, m)- 

To begin with let J7 C M" (n = 1, 2, 3) be a heat conductor with unknown cavities D C ft 
which consist of several disjoint domains such that D C ft and il\ D is connected. Each 
of these domains can be considered as a cavity. We assume that the boundaries dft and 
dD of n and D, repectively, are of class. Let the initial temperature distribution of 
n with cavities D be zero. By giving a heat flux / e ^^((0, T); H^^/^{dn)) on dfl over 
the time interval (0, T), we measure the corresponding temperature distribution on dft over 
(0, T) by an infrared light camera which is the case for thermography. If we can repeat this 
measurement so many times so that we can almost have the Neumann-to-Dirichlet map 

Ad : i'((0, T); H'^^^idn)) ^ ^^((0, T); H^^^{dn)) (1.1) 

which maps the heat flux / on dfl over the time interval (0, T) to the corresponding tem- 
perature distribution of dil over (0, T). This can be really the case if the given heat fluxes 
are confined at a part of dil and we have the cooling boundary condition at the other part 
of the boundary, because in this case the temperature generated by the heat flux goes to 
exponentially fast in time after we stop giving the heat flux. To make our argument simple, 
we will remain to the case there is no cooling boundary condition on the part of dfl. Then, 
our inverse boundary value problem is as follows. 

We define more precisely the Neumann-to-Dirichlet map Ad hy formulating a forward 
problem whose solution describes how the temperature distribution u in D,\D generated by 
the heat flux / diffuses. To begin with let 

W{{n\D)T) ■- [u e L^{{Q, T); H\Q.\D)), dtu G L^'dO, T); {H\n\D)y = H-\Q\D))j . 

Then, the forward problem is to look for a function u ~ — u{f ) G T4^((f2 \ D)t) which 

satisfies 

{{dt - A)u = in (Sl\D)-K (0, T) =: (ft \ D)t, 
d,u = fe L2((0, T); H-^/\dn)) on dfl x (0, T) {dft)T, 
d^u = OondDx (0, T) =: {dD)T, ' 
-u = at t = 0. 
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Here v is the unit normal vector to 9f2 (or dD) directed outside f2 (or U) and for any set 

E (1W,Et denotes the cyhndrical set x (0, T). 

This is a weU-posed problem, by which we can define the Neumann-to-Dirichlet map 

Kn : T); H-^^idQ)) ^ 1^(0, T); H^'\dn)), f ^ «/|(an).. (1.3) 

If there is no any cavity, i.e. = 0, we denote the Neumann-to-Dirichlet map by A0. 

When we take the Neumann-to-Dirichlet map Ajj as the measured data, our inverse 
boundary value problem is formulated as follows: 

Inverse Problem: Reconstruct D from A^. 

The linear sampling method considers some integral equation Fg = F^^ of the first 

kind for F =: Ajj — Ag, where F^^ is the Green function of the heat operator 9t — A in 
n X (0,r) with the Neumann boundary condition on dfl x (0,T) which has a singularity 
at {y,s). Then, the linear sampling method examines the behavior of an approximation 
solution to this integral equation for the case y G D and y ^ D when we fix s. A new feature 
which we can observe here is that there is a freedom to choose s, which suggests that we 
have more data than the linear sampling methods for the inverse boundary value problem 
for EIT and inverse scattering problem with near field data. 

The rest of this paper is organized as follows. In the next section, we will give the linear 
sampling method for the case of y G D and y ^ D, respectively. Some properties of the heat 
potential used in Section 2 will be proved in Appendix. 



2 Linear sampling method for identifying the cavity 

Let 

T^,,,^{x,t):=r{x,t;y,s)={ (M*-))-^ ("fe^) ' * > ^' 

I 0, t<s 

be the fundamental solution to the heat operator dt — A, and F°^ ^^(x, t) be the Green 
function of the heat operator with Neumann boundary condition on (90)t. We define the 
operators S, H, A, F as follows: 

• Sf := , Hf := d^S f\(^djj)^ with the solution to 

{dt - A)uf = in Ot, 
< dvU^ = / on {dfl)T, 
uf = at i = 0; 

• Ag := 2:^|(an)T ^i*^ the solution to 

{dt-A)zS = in {n\D)T, 
d„z^\(^OD)T = g, d^z^\(^oQ)T = 0> 
z'\t=o = 0; 
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• F ■.= Ad- Afl. 
Then, we have the following lemmas: 
Lemma 2.1 F = -AH. 

Proof. Let v = Sf with / e L^{{0, T); H~^/^{d^)) and consider the solutions w and u to 

{dt - A)w = in {n\~D)T, 
duw\(dD)T = -c>^v\(pD)^, = 0, 

w|t=o = 

and 

{dt - A)m = in {n\D)T, 

u\t=Q 0, 

respectively. Since {dt — A)u in (J7 \ D)t, we have 

(at - A)(w -f) = o in (17 \;d)t, 

_ {u-v)\t=o = 0, 

which implies w = u — v hy the uniqueness. Hence, by i9iyw|(aD)y = i?/, we have 

= {Ad - A^){d.vyn)^) = {Ad - A^)f = Ff. 
Therefore, F = -AH. □ 
Lemma 2.2 The operator 

H : l2((0, T); H-^'^{dn)) ^ L"((0, T); H-^'^{dD)) 
is continuous and has dense range. 



Proof. We prove the result by using layer potential argument for the heat equation. To 
this end, we introduce the function space H^''^ defined by 

H'-'^Tt) := i^((0, T); i?''(r)) n H'{{0, T); ^^(r)), T = or dfl. 



Recall that u := Sf satisfies 

{dt - A)m = in r^T, 
di^u — f on {dQ)T, 
u = at t = 0. 

We express the solution u to p.ip as a single-layer heat potential: 



(2.1) 



u{x, t) — KqC/) ■ 



Jdn 



r{x, t- y, s)(j){y, s)ds{y)ds, (x, t) e (2.2) 
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with an unknown density (f> G H'^'i ((c^ri)^)- Define the operator N by 

iV,/>:= i[7i(i^o0)k2+7i(^o0)|a=], 

where 

is a continuous linear map [5]. Then, using the jump relations of layer potentials in Theorem 
3.7 of [B], we derive that 

^.P + N<j>^f. (2.3) 
According to Proposition 4.3 in [Q, we know that 

i/ + TV : H^'i{{dn)T) ^ H^'^{{dn)T) (2.4) 
is an isomorphism. By defining 

Kt}} a^(i^oV')l(aD)T 
for V' e i2((o, T); R-^dn)), we have 

H = X(i + 7V)-i. 

In addition, the well-posedness of the forward problem says that the operator 

H : L^iO, T); R-^^^dn)) ^ L^iiO, T); H-^/\dD)) (2.5) 

is continuous. Note that ((951)^) is dense in L^((0, T); i?^^((9i7)). It is enough to 

show that the operator 

K : i2((o, T); H~^dn)) ^ 1^(0, T); H^^{dD)) 

has dense range. We give the proof as follows. 
By direct calculations, we have that 

Ki,{x,t)^ j I M{x,t;y, s)i^iy, s)ds{y)ds, {x,t)e{dD)T (2.6) 
Jo Jdn 

with 

M{x, t; y, s) := d„(^x)T{x, t; y, s). 

Let 

K* : L=^((0, T); H^dD)) ^ L^{{0, T); H^dQ)) 
be the adjoint oi K. It follows that 

K*viy, s)= [ I M{x, t; y, s)v{x, t)ds{x)dt, (y, s) e {dn)T. (2.7) 

JdD 
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For proving the denseness of K, it suffices to show that 77 = if 

{K^, ry) = (V', K*rj) = 

for aU iP e L^iiO, T); H-i{dn)). 
Define 

w{y, s)= [ [ M{x, t; y, s)rj{x, t)ds{x)dt, {y, s) G {M.'^\D)t. (2.8) 

Js JdD 

According to Appendix we have 

w in {M."\D)t, 

whence 

li; = on {dD)T- 
By the uniqueness result for the heat equation in D^, we also have 

w = in Dt- 

Consequently, we obtain that r/ = from the jump formula of the operator K* as y tends 
to dD, which is shown in Appendix |B] The proof is complete. □ 

We are now in a position to state the main result for the linear sampling method: 

Proposition 2.3 Fix s G (0, T) and let y Cz D. Then, there exists a function G 
L^((0, T); H^'^/^{dfl)) satisfying the inequality 

11-^5 - r^j;,s)llL2((o,T);ffi/2(an)) <e (2-9) 

such that 

^lim^ll5^llL2((o,T);if-i/2(an)) =00 (2.10) 

and 

J^^\\Hg^\\L^((o,Ty,H-^/HdD)) =oo- (2-11) 

Proof. We will adopt the proof given in [3] and PT|, which is for the scattering case. Let 
be the norm of the operator 

A : L^{{0, T); H-^'^{dD)) ^ L'^HQ, T); H^''^{dn)). 

For given e > 0, by Lemma [2.21 there exists a function g^ G -^^((0, T); H^'^/'^{d^)) such 
that 

^ 

WHg'" - (-c^i/r(^_^))||^2((o,T);H-i/2(aD)) < 
Since 9^r^^_ [(an)^ = 0, we have 



Hence, by Lemma [2. 11 we have 

\\Fgy - r°j^^5)||L2((o,T);//i/2(9o)) 
= II - AHgy - ^(9,.r^j^ ,,)|(9i3)^)||i2((o,T);_f/i/2(an)) 

< ll^ll 11^.9^ - i-di'^°y,s)\{dD)T)\\L^HO.T};H~^/^(dD}) < 

On the other hand, due to the boundedness of H, we have 

ll.9''llL2((0,T);_H"-i/2(ao)) 

> c||i?5^||^2((o^T);_ff-l/2(ai5)) 

> C (^||9^r°^^s)||L2((0,T):ff-i/2(a£))) ~ 11^.9^ - i-d^^"y,s))\\L^i(0,T);H~^/^(dD))j 

> 4d^'^°y,s)\\L^{{0,T);H-V2(0D)) - pjj 

with some constant c > 0. 
We next prove 

\\d,^'^°y,s)\\Lmo,T);H-^/2(dD)) ^ (y ^ dD). 

We only consider the case of n = 3 (for n = 1, 2, the proof is much easier). By the 
assumption on dD E C^, for any point xq G dD there exists a C^-function / such that 

D n B{xo, r) ^ {x e B{xo, r) : x^ > f{xi, X2)} , 

where B{xo, r) is the baU with radius r > and center at xq. We choose a new orthonormal 
basis {cj}, j — 1, ■ ■ ■ , 3, centered at xq with 63 — —v, where v is the normal unit vector to 
the boundary at xq. The vectors ei, 62 lie in the tangent plane to dD at xq. Let 77 be the 
local coordinates defined by the basis {cj}. Define the local transformation of coordinates 
1] = F{x) as follows: 

7]' = x', ri3 = X3- f{x') with x' = (xi, X2), rj' = (771, 772). 

Thus, in the new coordinates system, we can let ^ = (0, 0, ^3), ^3 > and if E [—1, I] x 
[—1, I] D2 with a sufficiently small constant I. Without loss of generality, we set s — 0. 
Then, we have 




where ci is a constant. 

We now estimate the last term in above inequahty. To this end, we introduce an auxiliary 
function (f{r]') = ce~~^ . By direct calculations, we obtain that 

Jd2 J's? 



Hence, we deduce that 



/ \Vr,'^{v')\'dr,'<c'^ f W\'e-'^dv' = c' [ IVpe-l-^'l^dy 
Jd2 Jm.'^ 

p poo pOO 

JR2 Jo Jo 



= 27rc2 
and therefore 



1 



= TTC 



Mm(D,) < {2y^t + tt) < {2y^T + tt) . 
In the following, we take c to satisfy 

(2V^r + 7r) < 1. 

Note that 

Hv) • - = (0, 0, -1) • (r?', -^3) = 6, 1^ - e = W\' + el- 



We have 



(^4^)3 2t 



exp 



4t 



At 



Taking (p = ce yields 



sup 

lvljfl(D2)<l 



J D2 



dt 



> 



167r3/2 
c 



D2 



87r3/2 



Jd, 



with 
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If we note that 

I{t) := / e-l^'I'dy 
is monotonically decreasing in [0, T], it can be concluded that 

ll'5i/r°^,s)lli2((o,T);/fi/2(a_D)) 



647r3 
c^cjlHT) 



^3'^ L2 se'^ds -^00 (^3 as 2/ dD). 



647r3 

The proof is complete. □ 

To further investigate the behavior of the density as y approaches to dD from the 
outside of £>, we need the following two lemmas. 

Lemma 2.4 The operator 

A : L^{{0, T); H-^'^{dD)) ^ L^HQ, T); H^''^{dn)) 

is injective and has dense range. 

Proof. The injectivity immediately follows from the unique continuation property for dt — A. 
For the denseness, we will adopt the proof of Lemma 4.1 given in [5], which is for the 
scattering case. Let gj e ^^^((0, T); H~'^l'^{dD)) {j e N) be such that the linear huU [{gj}] 
of {gj} is dense in £^((0, T); H^^/^{dD)). By Hahn-Banach's theorem, it is enough to 
prove the following: 

/ e i2((0, T);H-^''^{dn)) : [ ipjfdsdt = for all j G N with >pj := |(ao)T 

implies / = 0. 

Consider the solution v € W{{i} \ D)t) to 

{dt+A)v = in {n\D)T, 
v\t=T = 0, 



and set Zj — z^^ . Then, we have 



= 



/ [vAzj ~ ZjAv)dxdt 
J(n\~D)T 



'(0\D 

[dijZjV — dyVZj)dsdt — I [duZjV ~ d^vzj)dsdt 

{dn)T J{dD)T 

gjvdsdt. 

(dD)T 

Hence ti = on {dD)T- Using this fact together with = 0, we have w = in 

(57 \ D)t- Therefore / = 91^^1(911)^ = 0. The proof is complete. □ 
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Lemma 2.5 For any fixed s G (0, T), V'^^^ ^ R{A) range of A if y e Vt \ D . 

Proof. Suppose TO^^^ e R{A). Then there exist a function / e ^^((0, T); H-^/^{dD)) 
and the solution u;-'' e VK((r2 \ £')t) to 

{dt - A)wf in (r2\!D)T, 
w^t=o = 



such that 

Since 9i/W-'^|(an)T — '^•^^ly s)\{d^)T — Oj we have by the miique continuation principle 

w^^Tl^,) in (n\(:DUM))^ 

and hence 

\\^°y, s)\\ L^iiO.T); H^n\D)) = W^'^ W {{0,T); H^n\D)) < 

Here y £ rt\ D can be either y e 9-D or y e \ -D, but which ever the case it holds 

\\^°iv,s)\\L^{{0,T);Hi{n\D)) = oo- (2-12) 

In the following, we give a proof for (j2.12p . 

Let n — 3, s = and y e dD (The other case can be shown more easily). Then, for some 
S{0 < S < 7r/2), 7 (I7I < tt), i?, t (0 < i?, t < 1) and any e (0 < e < 1), we have 



647r3 



> 



^[-cos^]r'27^%-'dt^ r^exp (^-^) dr 
1 _ 

—7 cos (5 / t^^dt^ 00 (e ^ 0). 



Therefore 



and hence 



\\^ly,s)\\L^{{0,T);L^(n\D)) — °° 
\\^ty,s)\\L^({0,T);H^(n\D)) = OO- 

This is a contradiction. □ 



In contrast to Proposition [23] for y G I?, we now establish the following blowup property 
of the density for y ^ D as, follows: 
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Proposition 2.6 Fix s E (0, T) and let y £ D, \ D . Then, for every e > and 5 > Q, there 
exists g = .9" G -^^((0, T); H^^^'^{dfl)) satisfying the inequality 

W-r1y^4mio,Ty,m/^an))<e + S (2.13) 

such that 

1™ ||g^||L2((o,T);ff-V2(an)) =00 (2.14) 

and 

lim ||i?5^||^2((o,T);H-i/2(aD)) =00. (2.15) 

Proof. We will adopt the argument given in [3] and [T7], which is used for the scattering 
case. Since A is injective and has a dense range, for arbitrary S > there exists a function 

with the singular system (/j,„, ipn, g-n) of A such that 

II A/a - r°^,s)llL2((o,r);Hi/2(an)) < 5. 
By Lemma [2.51 we have T^^^^ ^ R{A). Hence, by Picard's theorem, we have 

ll/allL2((o,T);_fi--i/2(a_D)) ^00 (a ^ 0). 

By the denseness of R{H) in L2((o, T); H~^/'^{dD)), there exists 

g^eL2((0, r);i?-i/2(ai7)) 

such that 

ll^/a + ^^5allL2((o,T);_f/i/2(aa)) < £■ 

Hence, we have 

= II - ^Hg^ - r°^^s)||L2((0,T)://i/2(9o)) 

< II - AHgl, - ^/a||L2((0,T);//i/2(9o)) + \\Af^ - ,,)||i2((o^T);fl-i/2(an)) 

< £ + 5. 

To finish the proof, we recall that 

Hgl « fl in L2((o, T); 

and 

ll/allL2((o,T);_f/-i/2(9D)) OO (a ^ 0). 

These imply 

\\Hgl\\Lmo,Ty, H-^HdD)) -> oo (a ^ 0) 
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and hence 

ll5QllL2((o,T);ff-i/2(an)) oo (a -5> 0). 
Then, the rest of the proof can be done by noticing a = ol{8) — > ((5 — > 0). □ 

By Proposition 12.31 and Proposition 12.61 we have the hnear samphng method for identi- 
fying cavities in a heat conductor by measurements at its boundary. 
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A Property of potential (12.81) 



This Appendix is devoted to showing that 

(•T 



w(x, t) = / M{y, s; x, t)Tj{y, s)ds{y)ds = in (R" \ D) x (0, T) 

Jt JdD 

for It; = on {dft)T- It can be easily verified that 

dtw + Aw = Om (M" \D)x (0, T) 

and 

w = at i = T. 

Denote by N{K*) the nuUspace of the operator K*. Let i] = r/{x, t) € N{K*) be such 
that rj e C^{{dD)T) and ?7|t=o = 0. Note that such 77's are dense in N{K*). We continuously 
extend rj to t < and use the same notation for the extended one. 

We now establish the uniqueness of solutions to the following problem: 



dtW + Aw^O in (M" \ H) x (-00, T), 
w = on af2 X (-00, T), 
w = at < = T. 

For s > t and x ^ y, using Theorem 6.15 in |16) and the inequality 

r'^e-'^ < /3^e-^, < r, /3 < +00, 

we deduce that 



(A.l) 



\M{y, s; X, t)\ 



{v{y), x-y) 



< 



< 



2^47r(s - t) s - t 



■ exp 



(s - t)°'^\x - y|»-2ai 



■ exp 



\x-y\'^ 
'A{s-t) 



'4(s-t) 



(A.2) 
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for < ai < 1/2, where Ci, C2 are constants. Similarly, we also have 



du(x)M{y, s; x, t)\ 

1 {'^{y),x-y){iy{x),y-x) 



4y^47r(s - t) s - t 

1 iHy), ^(^)) 



s - t 



■ cxp 



4{s-t) 



< 



2^47r(s-t) s - t 

exp 



■ exp 



C, 



\x ~ y\' 



1 



\x-y\^ 

\x~y\^ 



exp 



< 



^ + 

(s - i)"2|a; _ y|"-2a2 (^s -t)°'-^\x -y\ 



4(s - t) 

\x~y\^ 
'4(s-i) 



2+n-2Q3 



(A.3) 



for < a2, 0:3 < 1/2, s > t and x ^ y, where Cj (j = 3, • • • , 6) are constants. 

Let Bfi be the disk with sufficiently large radius R and center at the origin. Define 
:= {W^\Ti)r\ Br. Then, by the estimate (|X2|) . we have 

-oo)IIl^(o«) = 0. 

Making use of this fact, combined with \\w{-, T)||i2(si^) = and w = on 9il x (—00, T), 
we obtain that 







Slnx(-oo,T) 



w{dt + A)wds{y)ds 



S1rx(-oo,T) 



wd^wds{y)ds. (A. 4) 



aSRX(-oo,T) 



By the estimates (jA.2p and (jA.3|) . we can easily deduce that 



/ 



It follows from (|A.4|) that 



(-00, T) 



lim 

R-¥oo 



wd^wds{y)ds — > (i? — > +00). 



\V.w\' = 0, 



Slj,x(-oo,T) 



and therefore 



w = in (R" \n)x (-00, T) 

by observing it; = on x (—00, T). 

Due to the unique continuation principle, we have 

u; = in {W \D) x (-00, T). 

The proof is complete. 
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B The jump relation of K* 



In this appendix, we show the jump relation of the adjoint of K which is used in proving 
Lemma 12.21 At first, we justify it for the continuous density by using the argument in the 
proof of Theorem 9.5 in [IS]. 

We interchange the variables (x, t) and {y, s) in (I2.7p . that is, 

K*7^{x, t)^ I [ M{y, s; x, t)7]{y, s)ds{y)ds (B.l) 

with 



t JdD 



Define 



w{x, t) := [K*vi){x, t) 

I Miy, s; X, t^y, s)dsiy)ds, {x, t) e (M" \ dD) x (0, T) (B.2) 

t JdD 



and set a = ttt^- Then, we have 



\x-y\ 
w{x, t) 



1 f {v(y), x-y) ( \y-x\ , n , / n , 

■^^P ~T? 7{]V{y-,s)ds{y)ds 



t [Ms-t)r/'JoD ^(s-t) "V 4(s-t) 

cr exp(-CT )77(2/, t H -—^)d(jds{y). 



JdD \x-yr J 4^ 4ct2 

Therefore, we can view the heat potential (|B.ip as a harmonic double-layer potential with 
density 



V-Ca:, y, := ^/ a""! exp(-a2)77(y, t + ' ^' 



4(72 

We now show that ip is continuous on M" x dD x [0, T) with 

lim ^p{x, y, t) ^ inViy, t) (B.3) 



for all y G 9D and < e [0, T), where 

7„ = / s"~^exp(-s^)ds. 
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To this end, we split tp{x, y, t) into three parts as follows: 



y, t) 



1 j^^/\x-v\ It-?;P 



2VT-t 



1 



r 1^ „,i2 



a"-^ exp(-a^) ^ 7y(y, t + ^-^) - 7y(y, t)) da 



1 Z"^"" 

+?7(y, / CT"-iexp(-CT2)ciCT 



:= h+h+h- 
Obviously, we have 

lim Ii{x, y, t) =0 
uniformly on dD and on compact subintervals of [0, T) and 

lim l3{x, y, t) = 7„77(2/, t) 

x-^y 

uniformly on dD x [0, T]. Since rj is uniformly continuous, for any e > there exists 5 > Q 
such that \ri{y, ti) - r]{y, ^2)! < £ for aU y e dD and all h, t2 G [0, T] with \ti - ^2! < S. 
Then, for all |x — ?/| < 45 and all a > y^\x — y\ we have 



viy,t+^-^^^)~viy, t) 



< e, 



and therefore 

|/2(x, y, i)| < / cr"~^exp(-CT2)dCT < 7„£. 

This means that 

lim l2{x, y, t) = 

x^y 

uniformly on dD x [0, T]. 

Using Theorem 6.13 in [16 , together with (IB. 31) . we can deduce that 

lim w{x =F hv{x), t) — ( [ M{y, s; x, t)7]{y, s)ds{y)ds ± —ri{x, t) (B.4) 
'«-*-o Jt JdD 2 

with uniform convergence on dD and on any compact subinterval of [0, T). This completes 
the proof for the case of continuous density. 

Based on the Lax theorem ([H]), the validity in can be shown by using the argument 
of Kersten in [2] . The proof is complete. 
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